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Abstract 

Relativistic mean-field models of nuclear structure have been enormously successful at repro- 
ducing ground-state properties of finite nuclei throughout the periodic table using a handful of 
accurately calibrated parameters. In this contribution we use powerful theoretical, experimental, 
and observational constraints — not employed in the calibration procedure — to validate two such 
models: NL3 and FSUGold. It is observed that FSUGold is consistent with all these constraints, 
except perhaps for a high density equation of state that appears mildly softer than required by 
astronomical observations. It is argued that incorporating such constrains goes a long way in 
removing much of the ambiguity left over from the standard calibrating procedure. 

PACS numbers: 21.65.+f,26.60.+c,21.30.Fe 
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I. INTRODUCTION 



Mean field descriptions of tlie ground-state properties of medium to lieavy nuclei liave 
enjoyed enormous success. Tliese fiiglily economical descriptions encode a great amount of 
physics in a handful of model parameters that are calibrated to a few ground-state properties 
of a representative set of medium to heavy nuclei. An example of such a successful paradigm 
is the relativistic NL3 parameter set of Lalazissis, Ring, and collaborators [H, [2|. 

Yet by their mere nature, such mean-field models are untested away from their narrow 
window of applicability. Whereas models fitted to the ground-state properties of finite 
nuclei all tend to agree on the saturation properties of symmetric nuclear matter, they 
widely disagree on its density and isospin dependence jl, |^. To resolve this ambiguity, an 
effort has been made to incorporate into the calibration procedure breathing-mode energies 
of heavy nuclei with neutron-proton asymmetries [b = {N — Z)/A] ranging from 6 = 0.11 to 
6 = 0.21 [5|,[6|. Although these approaches — combined with improved experimental data and 
analyzes — have helped narrow down the range of acceptable values of the incompressibility 
coefficient of symmetric nuclear matter to _ft' = 230±10 MeV 0, S, 0|, the density dependence 
of both symmetric-nuclear and pure-neutron matter remains largely undetermined. Thus, 
it is the aim of this contribution to test the validity of two accurately-calibrated relativistic 
mean-field models — NL3 [l|, [1] and FSUGold ^ — against recent theoretical, experimental, 
and observational constraints not employed in the calibration procedure. In this manner we 
aim to establish the extent to which models that were accurately calibrated around nuclear 
saturation density may be reliably extrapolated to the low- and high-density regimes. 

From the theoretical perspective, powerful arguments have provided critical insights into 
the behavior of pure neutron matter at low densities. Indeed, the low-density behavior of 
dilute Fermi gases with very large scattering lengths is universal, in that its energy equals 
that of the free Fermi gas up to a dimensionless universal constant of the order of 1/2 [l^, 11 



Yet pure neutron matter deviates from unitarity due to its relatively large effective range. 
Fortunately, effective-range corrections to unitarity were recently computed by Schwenk and 
Pethick [ij] . Such a model-independent approach will be used to test the validity of mean- 
field models away from their region of applicability. Moreover, it will be shown how such 
powerful theoretical constraints may be used to rule out a variety of accurately calibrated 
models 

Laboratory experiments with heavy ions have played a critical role in constraining the nu- 
clear equation of state. By tuning the energy of the colliding beams and the neutron-proton 
asymmetry, heavy-ion collisions probe vast regions of the phase diagram. For example, by 
compressing nuclear matter to pressures never before attained under laboratory conditions, 
the equation of state of symmetric nuclear matter was determined up to densities of about 
4-to-5 times that of normal nuclear matter 3]. Moreover, low-density constraints on the 



equation of state of neutron-rich matter are starting to emerge from the distribution of 
fragments in medium-energy collisions. In particular, by plotting the data in an ingenious 
manner, a powerful scaling relation — known as isoscahng [13, [l5| — was uncovered and 
shown to be sensitive to the low-density behavior of the symmetry energy [3. Note that 
the symmetry energy equals to an excellent approximation the difference between the energy 
of pure neutron matter and that of symmetric matter. 

Finally, enormous advances in both land- and spaced-based observatories have brought 
the fields of nuclear physics and astrophysics closer than ever before. In the particular case of 
neutron star structure and its intimate connection to the equation of state of dense matter, a 
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few recent developments are worth mentioning. Among these, observations of neutron-star- 
white-dwarf binaries with the Arecibo radio telescope have resulted in the largest neutron- 
star mass ever reported M(PSR J0751 + 1807) = 2.1 ±0.2 Mq [17|]. The limiting mass of 
a neutron star represents the optimal (and perhaps unique) way of constraining the high 
density component of the equation of state. If the above error bars can be narrowed down 
any further, a significant number of (soft) models of the equation of state will be ruled 
out 1^, 1^ . Moreover, transient phenomena — such as X-ray bursts powered by the nuclear 
burning of H/He in low-mass X-ray binaries — provide a powerful observational constraint 
on the long-sought mass-vs-radius relationship of a neutron star. For example, the mass- 
to-radius ratio of the neutron star in the low-mass X-ray binary EXO 0748-676 has been 
recently constrained by detecting the gravitational redshift of certain absorption lines 20 



More recently, by combining this result with additional observational constraints, both the 
mass and radius of the neutron star in EXO 0748-676 were individually determined to be 
M > 2.10±0.28 Mq and R > 13.8±1.80 km |2l| — an interpretation that appears to rule 
out all soft equations of state. For a comprehensive recent account on constraints on the 
equation of state from neutron-star observations see Ref. [18]. 

The manuscript has been organized as follows. In Sec. [TT] a brief review of the relativistic 
formalism will be provided. Particular emphasis will be placed on those model parameters 
that are loosely constrained by existent ground-state observables and the extent to which 
they may be firmly determined by the various constraints discussed above. Next, in Sec. IIIII 
results will be presented for a variety of low- and high-density observables that will be 
compared against available theoretical, experimental, and observational results. Finally, we 
offer a summary and conclusions in Sec. IIV[ 



II. FORMALISM 



The Lagrangian density employed in this work is rooted on the seminal work of Walecka, 




Serot, and their many collaborators (see Refs. 
first published by Walecka more than three decades ago [2 
implemented to improve the quantitative standing of the model. 



employ an interacting Lagrangian density of the following form j^, [IH], 2^ : 



and references therein). Since 
, several refinements have been 
In the present work we 



int 



The original Lagrangian density of Walecka consisted of an isodoublet nucleon field {ip) to- 
gether with neutral scalar (0) and vector {V^) fields coupled to the scalar density and 
conserved nucleon current [ip'j'^ip) , respectively jl^. In spite of its simplicity (the model 
contains only two dimensionless coupling constants), symmetric nuclear matter saturates 
even when the model was solved at the mean-field level \24!\. By adding additional contri- 
butions from a single isovector meson (6^) and the photon (A^), Horowitz and Serot 27 
obtained results for the ground-state properties of finite nuclei that rivaled some of the most 
sophisticated non-relativistic calculations of the time. However, whereas the two dimension- 
less parameters in the original Walecka model could be adjusted to reproduce the nuclear 
saturation point, the incompressibihty coefficient (now a prediction of the model) was too 
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Model 


TTT-s 


9l 


9l 




K 


A 


c 


Av 


NL3 


508.1940 


104.3871 


165.5854 


79.6000 


3.8599 


-0.0159 


0.0000 


0.0000 


FSUGold 


491.5000 


112.1996 


204.5469 


138.4701 


1.4203 


+0.0238 


0.0600 


0.0300 



TABLE I: Model parameters used in the calculations. The parameter k and the inverse scalar 
range rris are given in MeV. The nucleon, omega, and rho masses are kept fixed at M = 939 MeV, 
maj = 782.5 MeV, and mp = 763 MeV, respectively. 



large {K > 500 MeV) as compared with existing data on breathing- mode energies [28|] . To 
overcome this problem, Boguta and Bodmer introduced cubic {n) and quartic (A) scalar 
meson self-interactions that accounted for a significant softening of the equation of state 
(ir=150±50 MeV) [29]. Two parameters of the Lagrangian density of Eq. ([T]) remain to be 
discussed, namely, C, and Av. Both of these parameters are set to zero in the enormously suc- 
cessful NL3 model, suggesting that the experimental data used in the calibration procedure 
is insensitive to the physics encoded in these parameters. Indeed, Miiller and Serot found 
possible to build models with different values of C that reproduce the same observed proper- 
ties at normal nuclear densities, but which yield maximum neutron star masses that differ by 
almost one solar mass [25| . This result indicates that observations of massive neutron stars 
— rather than laboratory experiments — may provide the only meaningful constraint on the 
high-density component of the equation of state. Finally, the isoscalar-isovector coupling 
constant Av was added in Ref. |26| to modify the density dependence of the symmetry energy. 
It was found that models with different values of Av reproduce the same exact properties of 
symmetric nuclear matter, but yield vastly different values for the neutron skin thickness of 
heavy nuclei and for the radii of neutron stars [30] . The Parity Radius Experiment (PREX) 
at the Jefferson Laboratory promises to measure the skin thickness of ^osp]^ accurately and 
model independently via parity- violating electron scattering [sil, S^]- PREX will provide a 
unique experimental constraint on the density dependence of the symmetry energy due its 
strong correlation to the neutron skin of heavy nuclei [sj . 



III. RESULTS 

In this section we compare the predictions from NL3 [l|, 0] and FSUGold against the 
theoretical, experimental, and observational constraints discussed in the Introduction. The 
effective parameters of the two models are listed in Table [T] and their predictions for several 
bulk properties of nuclear matter are tabulated in Table HTl Note that po, £^o, and K denote 
the density, the binding energy per nucleon and incompressibility coefficient of symmetric 
nuclear matter while J and L denote the value and slope of the symmetry energy, all at 
saturation density. In particular, the pressure of pure neutron matter at saturation density 
is related to L through the following expression: 

^pnm(Po) = \poL . (2) 
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Model 


po(fm-3) 


eo(MeV) 


K(MeV) 


J(MeV) 


L(MeV) 


NL3 


0.148 


-16.24 


271 


37.3 


118.4 


FSUGold 


0.148 


-16.30 


230 


32.6 


60.5 



TABLE II: Bulk parameters characterizing the energy of symmetric nuclear matter (po, sq, and 
K) and the symmetry energy (J and L) at saturation density. 



A. Theoretical Constraints 

Perhaps surprisingly, one of the most stringent constraints on the equation of state of 
low density neutron-rich matter emerges from theoretical considerations, namely, from the 
universality of dilute Fermi gases with an "infinite" scattering length (a). In this limit the 
only energy scale in the problem is the Fermi energy [ep), so the energy per particle is 
constrained to be that of the free Fermi gas up to a dimensionless universal constant (^) 
that is independent of the details of the two-body interaction [13] • That is. 

To date, the best theoretical estimates place the value of the universal constant around 



While the neutron-neutron scattering length is large indeed (ann = — 18.5 fm), pure neu- 
tron matter deviates from unitarity due to a non-negligible value of the effective range of the 
neutron-neutron interaction (re = +2.7 fm). Thus, corrections to the low-density equation of 
state of pure neutron matter must be computed for kp ^rj ^ ~ 0.4 fm~^. Such corrections 
have been recently computed by Schwenk and Pethick [12], with their results displayed as 
the red region in Fig [TJ Also shown are the predictions of two microscopic models based 
on realistic two-body interactions, one of them being the venerated equation of state of 
Friedman and Pandharipande [s^. Finally, the predictions of NL3 and FSUGold are also 
shown. It is gratifying that the softening of the symmetry energy of FSUGold — caused by 
incorporating constraints from breathing-mode energies — appears consistent with the 
physics of resonant Fermi gases. Such a powerful universal constraint should be routinely 
and explicitly incorporated into future determinations of density functional. Indeed, such 
a constrain appears to rule out many of the models displayed in Fig. 2 of Ref. [§]. 



B. Experimental Constraints 

Laboratory experiments place important constraints on the equation of state of hadronic 
matter. Indeed, a variety of ground-state properties (primarily masses and charge radii) 
of finite nuclei are routinely incorporated into the calibration procedure of the models. 
However, the impact of heavy-ion experiments on these models is just starting to emerge. 
A particularly relevant example involves the distribution of fragments in medium-energy 
heavy ion collisions. It has been shown that the ratio of isotopic yields [R2i{N, Z)] obeys 



a powerful sca ling relation [IJ, Il5| that is sensitive to the low-density behavior of the sym- 
metry energy 16|. It is observed, quite naturally, that the reaction with neutron-rich nuclei 
{e.g., ^^^Sn+^^'^Sn) produces more neutron-rich and less proton-rich fragments relative to 
the neutron-deficient reaction {e.g., ^^^Sn+^^^Sn). This makes the reaction yields partic- 
ularly sensitive to the density dependence of the symmetry energy. For example, a stiff 
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FIG. 1: (color online) Equation of state of pure neutron matter as a function of the Fermi momen- 
tum. Predictions are shown for the accurately calibrated NL3 [l], 0] (green line) and FSUGold 0] 
(blue line) parameter sets. Shown also are various microscopic descriptions — including a model- 
independent result based on the physics of resonant Fermi gases by Schwenk and Pethick f]3] (red 
region) . 



symmetry energy (such as the one displayed by NL3 in Fig. [2]) imposes a stiff penalty on 
the system at high density for departing from the symmetric {N = Z) limit. It is, however, 
the softer symmetry energy (such as FSUGold in Fig. [2]) that imposes the stiffer penalty 
at the low densities relevant to the multifragmentation process. As such, one expects more 
neutron-rich fragments to be produced by a stiff rather than by a soft equation of state j36| . 
The experimental signature of this behavior is imprinted in a parameter than controls the 
variation of R2i{N, Z) with N for fixed Z (a parameter usually denoted by a). Using gen- 
eral thermodynamic arguments, the value of a was related to the change of the neutron 
chemical potential with the neutron-proton asymmetry 15| . Note that it was precisely 



the study of the neutron chemical potential of neutron-rich nuclei at zero temperature (i.e., 
the neutron Fermi energy) that lead to the conclusion that models with a soft symmetry 
energy reach the neutron-drip line before those with a stiffer symmetry energy 371 . 

In a recent study, Shetty, Yennello, and Souliotis [SS*] used the scaling behavior of the 
fragment yields — coupled to a molecular dynamics simulation — to extract the density 
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FIG. 2: (color online) Symmetry energy as a function of the baryon density (expressed in units of 
the saturation density po = 0.148 fm~^). Predictions are shown from the NL3 l2] (green line) 
and FSUGold 0] (blue line) models. Shown also are the results from Li and Chen 41] (maroon 
line) and the experimental analysis of Shetty and collaborators 3a] (red symbols). 



dependence of the symmetry energy. In the low density regime probed in the collisions, the 
density dependence of the symmetry energy may be parametrized according to the following 
simple formula: 

5(p) = 5(po)(^)' . (4) 

To make contact with this approach, our theoretical results were fitted to the above formula 
in the p = (0.3 — 1.0)po range, yielding values for the symmetry energy at saturation [5'(po)] 
and for the exponent (7) as displayed in Eq. ([5]). This same information is depicted in 
graphical form in Fig. [21 As in the case of pure neutron matter (see Fig. [T]) it appears 
that the density dependence of the symmetry energy predicted by FSUGold — at least at 
low densities — is consistent with this experimental analysis. However, the experimental 



determination is not without controversy [39|, |40| . Moreover, one should be aware that the 



connection between the collision of heavy ions and the zero-temperature equation of state 
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is model dependent. 



S{po 



( 31.6 MeV, Ref. [38], 

31.6 MeV, Ref. [4JJ, 

32.7 MeV, FSUGold, 
36.9 MeV, NL3; 
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'0.69, Ref. [38], 

0.69, Ref. [4JJ, 

0.64, FSUGold, 

^0.98, NL3. 



(5) 



"1 1 — 7] 1 1 r 

SNM 
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FIG. 3: (color online) Binding energy per nucleon as a function of baryon density (expressed in units 
of the saturation density on = 0.148 fm"'^) for symmetric nuclear matter. Theoretical predictions 
are shown for the NL3 [l|, 12] ( green line) and FSUGold 0] (blue line) models. Shown in the inset 
is a comparison between the equation of state extracted from energetic nuclear collisions fisl] and 
the predictions of these two models. 

Nuclear collisions may also be used to constrain the high-density behavior of nucleonic 
matter. To illustrate this point we display in Fig. [3] the binding energy per nucleon of 
symmetric nuclear matter as a function of the baryon density as predicted by both the 
NL3 and FSUGold models. Note that both models reproduce the equilibrium properties of 
symmetric nuclear matter and display the same quantitative behavior at densities below the 
saturation point. Yet their high-density predictions are significantly different. This emerges 
from a combination of two factors. First, FSUGold predicts an incompressibility coefficient 
K considerably lower than NL3, namely, 230 MeV vs 271 MeV (see Table [TTl)- Second, 
and more importantly, FSUGold includes a self-energy coupling [denoted by C in Eq. ([T])] 
that is responsible for a significant softening at high density. Note that the mixed isoscalar- 
isovector coupling (A^) plays no role in symmetric nuclear matter. We now compare the 
predictions of these two models against results obtained from energetic nuclear collisions 
that can compress baryonic matter to densities as high as those predicted to exist in the 
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core of neutron stars. The inset in Fig. [3] provides us with such a comparison. By analyzing 
the manner in which matter flows after the collision of two energetic gold nuclei, the equation 
of state of symmetric nuclear matter was extracted up to densities of 4-to-5 times saturation 
density [l3|. Figure [3] seems to rule out overly stiff equations of state (such as NL3). And 
while it continues to be gratifying that FSUGold is consistent with this analysis, one should 
reiterate that the connection between energetic nuclear collisions and the equation of state of 
cold nuclear matter is model dependent. Yet within these limitations, the same analysis [13] 
has been used to impose constraints on the equation of state of pure neutron matter by 
assuming two models — one soft and one stiff — for the unknown density dependence of 
the symmetry energy. The resulting equations of state are displayed in Fig. H] alongside the 
predictions of both theoretical models. The FSUGold parametrization appears consistent 
with a relatively soft symmetry energy. 
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FIG. 4: (color online) The pressure as a function of baryon density (expressed in units of the 
saturation density po = 0.148 fm~^) for pure neutron matter as extracted from energetic nuclear 
collisions IJ] by assuming soft (maroon) and stiff (red) symmetry energies. Also shown are the 
theoretical predictions from the NL3 [l|, y] (green line) and FSUGold 0] (blue line) models. 



C. Observational Constraints 

Ultimately, the most reliable constraints on the high-density component of the equation of 
state will come from astronomical observations. Indeed, important constraints are starting 
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to emerge from the combination of a large number of observations [18|. Here we limit 
ourselves to only two of them for their significant impact on the present analysis. For two 



other studies similar in spirit to the present one see the very recent references |42j and [43 

The first observation that impacts significantly on the present work is the one by Nice and 
collaborators at the Arecibo radio telescope [13]. Such observation of a neutron-star-white- 
dwarf binary system appears to suggest a neutron-star mass of M(PSR J0751 + 1807) = 
2.1 ±0.2 Mq (this is denoted by the red region in Figure E]). This appears to be the largest 
neutron-star mass ever reported and one that is significantly larger than those most ac- 
curately determined from double neutron-star binaries that display a mean of only about 
1.35-1.40 Mq with a very small dispersion. If the limits on the mass of PSR J075 1-1-1807 can 
be tighten any further (after all, at the 2a level the observation accommodates the rather 
wide range of 1.6-2.5 Mq) it would practically pin down the high-density component of the 
equation of state. Indeed, in the particular case of the FSUGold model with a prediction of 
only 1.72 Mq for the limiting mass, this observation would demand a mild hardening of the 
equation of state at high densities. As first observed by Miiller and Serot 2^ , in relativistic 
mean-field models this may be efficiently achieved by simply tuning the non-linear coupling 
while leaving intact all other parameters of the model. 




FIG. 5: (color online) Constraints on the mass-vs-radius relationship of neutron stars. Displayed 
in red is the allowed region as determined by the analysis of Nice and collaborators The black 
and orange solid lines result from the analyzes of EXO 0748-676 by Ozel 2l|], and Villarreal and 
Strohmayer 44 1 , respectively. Also shown are the theoretical predictions from the NL3 [l|, [3] (green 
line) and FSUGold jsl] (blue line) models. 



The second observation that seems to suggest a hard equation of state is that of the low- 
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mass X-ray binary EXO 0748-676. The first constraint on the equation of state from such an 
object came from the detection of gravitationally redshifted absorption hues in Oxygen and 
Iron by Cottam and collaborators [20|]. By measuring a gravitational redshift of z = 0.35, 
the mass-to-radius ratio of the neutron star gets fixed at M/i?~0.15 (with M expressed 
in solar masses and R in kilometers). By incorporating additional constraints arising from 
Eddington and thermal fluxes, a recent analysis by Ozel seems to place simultaneous limits 
on the mass and radius of the neutron star in EXO 0748-676. That is, M>2.10±0.28 Mq 



and R > 13.8 ±1.80 km 2l|. These limits are indicated by the black solid line in Fig. [51 
An earlier determination of the spin frequency of the same neutron star by Villarreal and 
Strohmayer [i^], when combined with the rotational broadening of surface spectral lines, 
yields an independent determination of the stellar radius of R^ll.b^H This estimate, 
when combined with the gravitational redshift, yields the orange line in Fig. [51 Note, 
however, that the use of rotational broadening to constrain the stellar radius has been put 
into question in Ref. j2l|. Finally, mass-vs-radius predictions from the NL3 and FSUGold 
models are displayed in Fig. [3 The results clearly indicate the significantly harder character 
of the equation of state predicted by NL3 relative to FSUGold. This, even when both models 
predict practically identical properties for existent ground-state observables of finite nuclei. 
Based solely on these observations, NL3 with its stiff equation of state appears to fair far 
better than FSUGold, despite the fact that both breathing-mode energies and heavy-ion 
experiments seem to suggest that the NL3 equation of state is overly stiff. 



IV. CONCLUSIONS 

Accurately calibrated relativistic models of nuclear structure have been enormously suc- 
cessful at describing a variety of ground-state properties throughout the periodic table by 
employing a relatively small number of effective parameters. Chief among these is the NL3 
parameter set of Lalazissis, Ring, and collaborators [1, 2]. Yet under closer scrutiny, it was 
revealed, perhaps not surprisingly, that the properties employed in the calibration procedure 
of such models are insufficient to firmly pin down the equation of state even around satu- 
ration density. In an effort to lift this "degeneracy" the FSUGold model was conceived j^. 
Relative to NL3, FSUGold includes two additional effective parameters. One of them [de- 
noted by Av in Eq. ([T|)] is responsible for a softening of the symmetry energy and should be 
firmly constrained by the pioneering Parity Radius Experiment at the Jefferson Laboratory 
that is scheduled to run in the early part of 2009 [32]. The other parameter [denoted by ^ 
in Eq. (Q] controls the high-density component of the equation of state and may be tuned 
to generate maximum neutron star masses that differ by almost one solar mass without sig- 
nificantly affecting the behavior of the equation of state near saturation density j25| . These 
two parameters — that are set to zero in the NL3 model — have a particularly dramatic im- 
pact on two observables: (a) the neutron-skin thickness of ^°^Pb and (b) the limiting mass 
of a neutron star. In both cases NL3 predicts significantly larger values than FSUGold, 
namely, 0.28 fm vs 0.21 fm for the former and 2.78 Mq vs 1.72 Mq for the latter. The main 
goal of this contribution was to test the validity of these models away from their region of 
applicability. 

The validation of FSUGold — a model calibrated to various ground-state properties and 
collective excitations of finite nuclei — was implemented through a detailed comparison 
against recent theoretical, experimental, and observational constraints. These constraints 



emerged from the universal behavior of dilute Fermi gases with large scattering lengths [12 
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heavy-ion ex per iments that probe both the low- and high-density domain of the equation 



of state [13|, [ij, ll5|, |38| . and astronomical observations that place limits on masses and 



radii of neutron stars [17|, l21|. On the basis of these comparisons, it was concluded that 
FSUGold meets all the challenges, even when no attempt was ever made to incorporate these 
constraints into the calibration procedure. If at all, only the observational data seems to call 
FSUGold into question by suggesting a slightly harder equation of state. The promise of 
improved observations and analyzes with existent and future missions, such as Constellation 
X, offers the greatest hope for determining the high-density component of the equation of 
state. 

The response to such future achievements should not be limited to an indiscriminate ad- 
justment of parameters. After all, "knobs" can always be turned to reproduce a particular 
set of experiments and/or observations. Rather, models should aim at reproducing simulta- 
neously a myriad of observables that probe the equation of state over a wide dynamic range. 
This should be one of the primary missions of all modern theoretical approaches. Only by 
adopting such strict standards one could test the validity and applicability of such models. 
Only then can one ensure the emergence of exotic phenomena. 
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